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ABSTRACT 


This method includes detecting the passage before a sensor 
of each of the teeth of the inertial flywheel of an engine; 
measuring the time of passage d, of each of these teeth; 
constituting a measurement horizon extending on both sides 
of the angular range of the explosions; dividing this mea- 
surement horizon into eight sets of teeth of equal length of 
which six are for the angular range of the explosions; 
measuring the passage times of these eight sets and in 
assigning them a rank in the measurement horizon du +1 , duo, 
. . . , du 5 , du +1 ; calculating the period of the explosions 
Todu 0 +du +1 +du 2 +du 3 +du 4 +du 3 ; calculating a term 
Q c o2du 0 -du_ 1 -du 3 -du +1 +2du 5 ; calculating a term D c op Q c 
wbere p=^/2V3, a constant stored in a memory; and calcu- 
lating the desired torque using the relation Cg c =A D c /T 3 + 
B/T 2 in which A and B are constants stored in a memory. 

20 Claims, 5 Drawing Sheets 
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INTERNAL COMBUSTION ENGINE described in the patent referenced above, or of the method 

TORQUE MEASUREMENT DEVICE AND and device according to the first object of the present 

METHOD invention, to measure the variations at a low frequency of 

the resistant torque applied to the engine to allow the 

BACKGROUND OF THE INVENTION 5 assignment of a quality factor to the torque measurement 

performed. 

1. Field of the Invention A third ob j ect of the mver)t i on concerns a method and 
The invention concerns a method and a device for raea- device derived from the preceding ones to correct the 

suring the torque of an internal combustion engine and, more low-frequency perturbations undergone by the engine 

precisely, to improvements to this method and device for 1Q torque, especially those caused by the oscillations of the 

measuring the aforementioned torque as described in French vehicle or by sudden changes in load. 

Patent No. 9,111,273 filed by the applicant. A fourth object of the invention is to propose a method 

2. Discussion of the Background and device to carr y out a new correction integrating the 
r ,. , . . , , . . . , effects of the richness of the mixture of the recirculation rate, 
In this patent a device is described to produce a numerical Qf ^ exhaust and of ^ transienl ■ conditions 

value Cg which is representative of the mean gas torque is q{ ^ mg ^ Qew cmelioa or DOt be 
generated by each combustion of the gaseous mixture in the combincd ^ the above-mentioned corrections, 
coders of a combustion engine at low level operating Accordiog tQ & ^ ^ q{ ^ inventi & method fof 
conditions. This device comprises: me calculation in a simplified manner of the term D men- 
measurement references arranged on a ring gear which is tioned above consists in- 
integrally connected to the inertial flywheel of the 20 cstablishing a measurement horizon equal t0 at least the 
engine or the shaft; angular range between consecutive explosions in the 
means to define at least one reference for indexing the engine; 

reference; regrouping the instantaneous passage time di of the ref- 

a sensor which senses the passage of the references and 25 erences contained in this measurement horizon in a 

which is mounted in a fixed manner in the vicinity of relatively small number of du 0 , . . . , du„ passage times 

the ring gear; of reference sets and in assigning to these times du their 

calculation means to produce a primary numerical value rank 0 ... m in each measurement horizon; 

d„ which is representative of the instantaneous passage combining together, by addition and subtraction, possibly 

time before the sensor of each of the references; 30 weighted, a number of times having predetermined 

calculation means to derive from the primary numerical ranks, so as to produce a magnitude Q which has a 

values d t a secondary numerical value T which is nonzero frequency response at the analysis frequency 

representative of the total passage time before the of the explosions of the engine; 

sensor of each series of n references in the angular determining a weighting constant such that the frequency 

range of the combustions of the engine; 35 response of the term D«p.Q is at least for said analysis 

calculation means to derive a second secondary numerical frequency and for the combination of the measurements 

value D, which is representative of the projection onto ( for Q. ,den " cal . * obtamed for the term 

a phase reference line of the references, corresponding calculated by said initial relation, 

to the origin of the angular periods of the combustions, According to a first embodiment of the preceding method, 

of the amplitude of the alternating component of the 40 the measurement horizon exactly corresponds to the angular 

instantaneous passage limes di of the references before interval separating two consecutive explosions of the 

the sensor, at the frequency of the combustions in the ""g^ the number of refereace bein g six and the six 

engine. The defining relation D is: measurement times referenced du 0 , . . . , dus, the magnitude 

0 Q=(du 0 +du s -du 2 -du 3 ) and the constant p=Jt/2V3, when the 

45 reference sets have the same angular length n/6. 

D - I di- cos(/ ■ 2*/n) Thanks to this simplified calculation method, the value 

obtained for the term D is essentially identical to that 
initially obtained, and this applies in spite of the fact that the 
hereafter called the initial relation number of arithmetic 0 p era ti 0 Qs to be performed has been 
calculation means to derive the desired numerical value 50 reduced to four. This considerably reduces the specifications 
Cg from the relation Cg*»A-d d/T^+B/T 2 , where A and 0 f the microprocessor used to perform the calculations. 
B are experimentally determined constants. Moreover, it should be immediately noted that the practical 
In a variant, the same patent proposes, to decrease the embodiment forms of the methods according to the inven- 
calculation times, to regroup the n references of an angular tion are determined by a systematic analysis of the math- 
range of combustions into three or four sets and to calculate 55 craatical equivalents of the initial relation which defines D. 
the above-mentioned term D from the passage times of these i n a second embodiment of the method, the measurement 
sets before the sensor. horizon is equal to the above defined horizon, increased by 
SUMMARY OF THF INTVPNTTDNT 30 additional P ackel du 6' lhe term °- **ing e <l ual 10 

SUMMARY OF THE INVENTION Q«du 0 -du 3 43/2. (du 6 -du 4 ). The second form has the advan- 

A first object of the invention is to propose a method and eo tage of improved discrimination between late combustions 

device for simplified calculation of the term D, using a small and combustion misses of the engine, 

number of reference sets and producing a result which is According to the second object of the invention, a method 

essentially identical to thai supplied by the initial calculation to measure the variations 8Cg of the resistant torque applied 

procedure, while at the same time considerably decreasing to a two- or four-cylinder engine consists in: 

the number of mathematical operations required. 65 establishing a measurement horizon which is equal to at 

A second object of the invention is to propose a method least the angular range between two consecutive explo- 

and complementary device either of the fundamental method sions in the engine; 
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regrouping the instantaneous passage times di of the Cg and that the terms k and J above are connected, the 

references contained in this measurement horizon in a variation of the load torque can be written, with good 

relatively small number of diio, . . . , du m , passage times precision: 

of reference set and identify these times du by their 6 Cg .^ 2 V3y(du 0+ du 5 -du ^du^AO*. 

rank 0 . . . m' in each measurement horizon during the 5 

course of an analysis, the angular values of the time sets Any nonzero value of this term oCg is a factor which 

of ranks 0 and m' being equal; negatively affects the quality of any measurement of Cg 

assigning the reference du.! to the time of the last performed at that instant. . 

reference sets of the last measurement horizon analyzed F ™ m tbe vahie of b ^ * * jn addition possible to directly 

before and the reference du +1 to the time of the first io apply a correction to the term D indicated above and to write 

reference sets of the next measurement horizon to be (a=p=JV2V3): 

analyzed; D-P(Q+q)*^vy(2du 0 -du_ 1 +2du s -du 3 -<tu 1hl ) 

combining four of these data elements du using the and 

relation q=flu=(du 0 -du_ 1 +du m .-du +1 ); 

deriving the desired term 6Cg using the relation 5Cgo 15 c ge *ADjr+Bn*. 

a.q.A/T 3 , where a is a new constant which is dependent Under these conditions, by means of a simple comple- 

on the angular value of the four reference sets con- mentary calculation, the corrected mean gas torque Cg. 

cemed. calculated from tbe corrected term 0 C is insensitive to 

According to a first embodiment of this method, the variations at low pressure of the resistant torque for as long 

angular length of each reference set is ji/6 and the constant 20 as the approximation of this torque by a parabola remains 

a=7t/2V3. valid within the measurement horizon. This is the case for as 

According to the third embodiment of the invention, a long as the value of the term 6Cg remains relatively low, 

method to produce a value of the term D defined above, with which corresponds to the absence of either large perturba- 

corrections for low-frequency perturbations which affect the tions of the resistant torque applied to the engine or abrupt 

engine because of changes in resistant torque applied to this 25 changes of large amplitude affecting this torque. The term 

engine, consists in calculating this value according to the 8Cg can quite clearly be produced independently of the term 

relation D. In this case, it constitutes a negative quality factor for the 

Dc=pQ+aq, where a and p are two constants which may torque Cg, regardless of the manner in which this torque is 

be identical. 3D calculated. 

The implementation of this third method is particularly According to the fourth object of the invention, the 

simple and effective, due to the considerations indicated method used to take into account tbe particular conditions of 

below, a combustion affecting the determination of the numerical 

The calculation of the mean gas torque Cg expressed by value Cg of the mean gas torque generated by each com- 

the fundamental relation Cg=A-D/T 3 +B/T 2 remains valid for bustion of gaseous mixture in the cylinders of a four-stroke 

as long as the load applied to the engine is decoupled for the and four-cylinder combustion engine at low-level operating 

analysis frequency considered. However, at low levels of conditions, and consequently to correct this numerical value, 

rotation, the analysis frequency approaches the oscillation which is expressed by the relation: Cg-AD/r 3 +B/T*, is 

frequency of the load. In tbe case of strong variations in the characterized in that it consists in: 

load (sudden demand for torque, oscillation of the vehicle on 4Q determining tbe air mass M pumped in a cylinder during 

a road in poor condition, initiation of transmission the period of combustion of another cylinder; 

resonance), the calculation of the mean gas torque is per- determining the pumping torque Cp required for this 

turbed according to the relation: purpose; 

Corrected Cg-measured Cg -kdQ/de, where C c is the storing in memory these two data element during the 

resistant torque applied to the engine (the load) and k, a 4J following two combustion periods; 

constant equal to the ratio of the mean gas torque to the combining this data to produce a correction term H„., for 

component of the instantaneous gas torque at tbe analysis the torque Cg ^ wm ch p Crta ins to the combustion 

frequency. In the case where the measurement horizon is of rank „. accor d m g to tbe relation: 
one-half turn or a complete turn, the limits of this horizon 

correspond to dead center for all cylinders. The torque 5Q ",'-M M ,-A-J 4flC / ,, '.2i 
exerted on the shaft by the pressure existing in the cylinders 

is thus zero since the lever arms are zero at that instant. The in which 8 and a are constants which depend on the 

load applied to the shaft can then be calculated by measuring type of engine, and 

the acceleration of the flywheel about dead center. calculating the corrected mean gas torque, generated 

Let d-and d+be the times corresponding to an angle d r on 5J during the combustion period of rank n', according to 

both sides of dead center. The load applied to the flywheel the relation Cg^.-bCg^+H,,., where h is a calibration 

is then: C c oJ7c 3 /de 2 (d+-d-)Ar 3 , where J is the total inertia constant which is dependent on the type of engine, 

of the flywheel and the shaft. By using C c _ and C e + to According to a particular characteristic of the method 

designate the torques so calculated respectively at the begin- according to tbe invention: 

ning and tbe end at the end of a measurement horizon, the 6Q determination of the air mass M consists in measuring the 

instantaneous variation of the load torque affecting the pressure Pc in the air-intake manifold of the engine and 

measured torque Cg is in calculating M-Pc-Rjv, where R„ is the essentially 

constant filling coefficient of the cylinders of the 

Corrected Cg-measured Cg-kip^-C^)!*. 65 determination of the pumping torque Cp consists in 

From the times of the sets with angular length ji/6 and, measuring the atmospheric pressure Pa and in calcu- 

taking into account that the term A of the relation defining lating Cp=(Pc-Pa). 
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The implementation of this particular characteristic of the 
method according to the invention makes it possible to 
express the term H„. according to the relation: 

Due to these arrangements, a correction is applied to the 
measurement of the torque of a four-stroke engine with four 
cylinders, which correction is adapted to the operation with 
a lean mixture, to a high recirculation rate of the exhaust 
gases and/or to the transitory conditions of the engine 
(variation in the butterfly valve angle). All these advantages 
result from the considerations indicated below. 

The torque exerted by the pressure of the gases on the 
crankshaft is the resultant of the torque applied by the 
cylinder in the expansion phase and of the torque applied by 
the cylinder in the compression phase. In the case of a 
theoretical thermodynamic cycle, the two expansion and 
compression torques are adiabatic, and thus identical to the 
nearest factor k (k increases when the combustion energy, 
and thus the richness, increases). 

In the case of a four-stroke four-cylinder engine, if C comp 
(6) represents the instantaneous torque exerted by the cyl- 
inder in the compression phase, the torque exerted by the 
cylinder in the expansion phase is then equal to: 

c exp-A-C^K-U). 

The gas torque is the sum of the expansion torque, the 
compression torque (negative because resistant) and of the 
pumping torque. Pc being the pressure in the intake manifold 
and Pa the atmospheric pressure, it is possible to consider 
the pumping torque to be proportional to (Pc-Pa). The gas 
torque is thus: 

C t „^-C^ p 4aiPc-Poy(k-l>C CMP +a(c-P a ). 

In the case of a computer estimation of the torque Cg 
defined by the relation indicated above, the symmetry of the 
curve of the coefficients with respect to the median plane of 
the measurement horizon is transformed into a curve with 
odd symmetry with respect to the torque. The result is that 
the resistant compression torque is considered a positive 
torque, just like the expansion torque, by the method for 
measuring Cg. The measured value of the torque is: 

Cg-ADfP+BrP-ViC^+c^. 

However, the compression torque C comp can be calculated 
from the air mass introduced into the cylinder, which can be 
derived, either from a measurement with a flowmeter or 
from a compression measurement. In the case of the 
pressure, the compression torque is proportional to the 
pressure Pc in the manifold and to the filling coefficient R N 
of the cylinder of the engine (which depends on the oper- 
ating condition of the engine): 

In the case of the generated gas torque, the torque C comp 
is the compression torque of the cylinder which will be in the 
expansion phase at the time of the measurement (pressure 
Pc c „_2), whereas for Cg, the measured gas torque, the torque 
C comp of the cylinder in the compression phase is measured 
at the time of the measurement (pressure Pc,„_ a ). In contrast, 
the filling coefficient can be the same for the two cylinders 
because there is little variation in the operations condition of 
the engine. By combining the above equations, we get: 
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or 

Cg c -h Cg+H, where h is a calibration constant which is a 
function the type of engine, and: 

5 H-bR N - (Pc^+Pc^+alPc^-Pa), 

in which a-V/jt, where V is the unitary displacement of 
the engine, and 8= 1.14V in a first approximation (with 
k=4.5) adjustable for each type of engine and filling 
10 coefficient R N filling coefficient of the engine, which is 
dependent on the operating conditions but is equal to 1 
in a first approximation (it varies from 0.8 to 1.02 as a 
function of the increasing operating conditions). 
It is recalled that the term R^.PC represents the air mass 
15 introduced into the cylinder. Consequently it can be replaced 
by another Da/N estimator, where Da is the flow rate of the 
air mass in the manifold and N is the number of rums per 
second of the crankshaft. Similarly, the term (Pc-Pa) can be 
replaced by (DT-Mq) , where M 0 is the air mass introduced 
2D into a cylinder for a fully-opened butterfly valve. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The characteristics and the advantages of the invention 
will become more precise in the description of a particular 
25 embodiment of the invention given below as a nonlimiting 
example, with reference to the drawings in the appendix, in 
which: 

FIG. 1 is a diagram of the different elements which 
constitute the measurement device according to the inven- 
30 tion; 

FIG. 2 represents the frequency responses of the two D 
terms, calculated by means of the initial relation and the 
simplified relation, respectively; 
35 FIG. 3 represents four implementation and operation 
waveforms of teeth-references packet times; 

FIG. 4 represents the diagram of the different elements 
which constitute a second embodiment of a device for 
measuring the torque of an engine, according to the inven- 
4Q tion; and 

FIG. 5 represents the diagram of the elements specific to 
a third embodiment of the invention. 

DESCRIPTION OF THE PREFERRED 
45 EMBODIMENTS 

According to FIG. 1, a circuit 10 is represented, for the 
measurement of the mean gas torque Cg^ corrected for 
low-frequency perturbations, which is produced by each 
combustion of the gaseous mixture in a four-stroke four- 
50 cylinder combustion engine equipped with a toothed mea- 
surement ring gear 12 which is integrally connected to the 
flywheel of an engine with electronic ignition. For example, 
the ring gear 12 includes, a its perimeter, fifty-six identical 
teeth which are regularly spaced, such as those formed by 
55 the solid part 14 and the recessed part 16, distributed into 
two series of twenty-eight teeth separated by two diametri- 
cally opposite reference teeth, such as the ring gear formed 
by the solid part 18 and the recessed part 20 which have a 
width which is twice that of the other teeth. In fact, the ring 
60 gear 12 comprises (2x28+2x2)=60 equidistant references 
which are made up of real or virtual teeth of the same 
module. The two broad teeth which are diametrically oppo- 
site are used as reference or origin of indexation to allow the 
numbering of each of the teeth and especially the identifi- 
es cation of the tooth do which will be defined below. 

With the four-stroke four-cylinder combustion engine 
referred to in the example above, it should be noted imme- 
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dialely that the angular period of the combustions concerns following stages the six reference times du 0 , .... du 5 of the 

30 teeth and it js equal to half the period of rotation of the central series as well as, under the reference du_j, the time 

crankshaft. du 5 of the last set of the preceding series and, on the 

A fixed sensor 22 is associated with crown 12, for reference du +J , the time du 0 of the first set of the series 

example, one with variable reluctance, which is adapted for s which has just arrived. This is represented in line II of FIG. 

the delivery of an alternate signal 24 having a frequency 3 - 

proportional to the speed of passage of the teeth of the ring A first calculation step 42 which receives these eight 

gear that is proportional to the instantaneous speed of the successive times is programmed so as to select the six 

flywheel. reference times du 0 , . . . du 5 and to calculate from them the 

The angular position of this sensor 22 with respect to the 10 sum, which is the instantaneous period T of an angular range 

indexing teeth 18 at the time when the piston of a cylinder of combustions. Each new value T so calculated is addressed 

is at top dead center is known or determined. This allows the to « buflfcr memory 44 where it replaces the value calculated 

identification of the tooth do as being the one which passes earlier. 

before the sensor during the passage of the piston of a given A second calculation step 46, which receives these eight 

cylinder through its top dead center of combustion. The 15 successive times (du^, du<j, . . . , du 5 , du +1 ) is programmed 

signal delivered by the sensor 22 is applied to the inlet of a to select the four following times: d^, di^, du 3 and duj, and 

shaping circuit 26 which is adapted to the delivery of signals du 5 , and to combine them according to the relation repre- 

28 with steep edges, having a period equal to the instanta- sented with the signs and the weightings in line II of FIG. 3: 
neous period d,- of the incident signals 24, the index i varying 

from 0 to 29 as the teeth pass before the sensor. Each period 20 Q-^-^au^). 

d, of the signal so produced corresponds to the passage time a third calculation stage 47 which receives the same eight 

of one tooth, either a solid part or a recessed part, before the times is programmed to select the four time du_ 1? du 0 , du 5 

sensor 22. As far as the incident signals produced by the and du +1 and to combine them according to the relation: 

reference teeth 18 are concerned, the shaping circuit 26 q=6u=(du 0 -du_ 1 +du 5 -du +1 ). 

transforms them in the same manner into a signal with steep 25 The terms Q and q so calculated are addressed to a buffer 

edges, having a time which is exactly twice that of the memory 48 where they remain until new values of Q and q 

signals pertaining to the other teeth. The signals 28 are arrive. The buffer memory 48 is connected to a fourth 

applied to a stage 30 for the measurement and the calculation calculation stage 50 which also receives from a memory 52 

of the instantaneous periods d, of passage of the real and the weighting constant p=a=rt/2V3. The calculation stage 50 

virtual teeth of the measurement ring gear 12 before the is adapted for calculating the sum Q c -Q+q, and then the 

sensor 22. product D c =p-Q c , and to apply it to a fifth calculation stage 

The measurement and calculation stage 30 comprises 54. 

counting circuits which receive high-frequency timing This calculation stage 54 also receives from a memory 56, 

pulses (10 MHz, for example) produced by a quartz clock 32 35 two constants, A and B, experimentally determined and, 

and it delivers at the output numerical values which are from the buffer memory 44, the term T. From these four 

representative of the number of clock pulses counted magnitudes, me stage 54 calculates the mean gas torque, 

between two low-high transitions of the signals with steep corrected for low-frequency perturbations according to the 

edges produced by the shaping stage 26. As far as the relation Cg c =AD c /T 3 +BAT 2 . As a consequence of what has 

processing of each signal 28 produced by one of the long 4Q been said above, it should be noted that this value so 

indexing teeth 18-20 is concerned, in the context of the calculated applies only to two- or four-cylinder engines, in 

described example, its value will be divided by two (one-bit which high or low dead center occurs with every turn or 

shift) , and the result will be assigned to the corresponding every half turn of the crankshaft, depending on whether the 

virtual teeth. In this manner, the stage 30 for the measure- engine is a two-stroke or four-stroke engine. For engines 

ment and calculation of the periods d,- sends a successive 45 with six or eight cylinders, for example, the correction of the 

series of thirty numerical values d t - to a buffer memory 34, mean gas torque for at low-frequency perturbations will no 

which are respectively associated with thirty consecutive longer be possible. In this case, the buffer memory 38 will 

numbers ranging from 0 to 29 defining the rank i of each of be limited to a series of six sets of teeth-references and the 

the real or virtual teeth within a measurement horizon. calculation stage 42 will be adapted so as to perform only the 

The buffer memory 34 is connected to a calculation stage 50 calculation of CKdu 0 -du 2 -du 3 +du s ). 

36 which is adapted for regrouping the thirty instantaneous With Q c -Q+q, the correct mean gas torque Cg^, calcu- 

times d,. between two reference teeth in sets of fine. In this l ated ™ Seated above, comprises two elements, that is the 

manner, the stage 36 successively produces six passage torque Cg measured according to the simplified relation 

times of five teeth -reference sets, respectively bearing the usin 8 the term 

numbers di^, . . . , du 3 as they are produced. These six times 55 Q=(du 0 -du2-du 3 4.du 3 ) and a correction 5Cg applied to 

and their six rank numbers are represented in line I of FIG. this torque Cg so measured, using the term 
3. These times and these numbers are then transmitted to a 

buffer memory 38 which is placed under the control of a ^o-^-i-^i+^s)- 

transfer control stage 40 which is also adapted for receiving Consequently, the precision of the measurement of the 

the times du 0 , . . . , du 5 , and identifying the rank of the time 60 torque Cg, which is carried out from the simplified relation 

du 0 . using the term Q, is replicated entirely in the measurement 

The buffer memory 38 is adapted to contain three sue- of the torque Cg. corrected for low-frequency perturbations, 

cessive series of six times of sets-references: a central series In contrast, this is not the case with perturbations at fre- 

which defines the measurement horizon during the analysis, quencies which are much higher than the frequencies of the 

a front series which precedes the central series in time and 65 combustions. For these high-frequency perturbations, nei- 

a rear series which follows it. When a new time du 0 arrives, ther the initial calculation procedure nor the simplified 

the transfer control stage 40 is adapted so as to transfer to the calculation procedure according to the first object of the 
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present invention allows their correction, the simplified 
calculation procedure being in that case worse than the 
initial calculation procedure. These different results clearly 
appear in FIG. 2 which represents the frequency response of 
the two terms D calculated respectively according to the two 5 
procedures in question. Indeed, these frequency responses 
are identical up to four times the frequency of a measure- 
ment horizon 2Qo (Qo being the frequency of rotation of a 
four-stroke, four-cylinder engine). At higher frequencies, the 
frequency response of the term D, obtained using the sim- io 
plified calculation procedure, presents amplitude maxima 
which are much greater than those of the term D obtained by 
the initial calculation procedure. This signifies that the 
high-frequency perturbations are taken more into account by 
the simplified calculation procedure than by the initial is 
procedure. 

According to this FIG. 2, it is verified that the magnitude 
Q has a nonzero frequent frequency at the analysis fre- 
quency of the explosions of the engine, and a nonzero 
frequency response at zero frequency, that is a zero value for 20 
a strictly constant speed of rotation of the engine, which is 
essential for the quality of the approximation D=pQ. 

In an embodiment variant of the method, the measure- 
ment horizon is equal to the horizon defined earlier, 
increased by an additional set du e , for example, du+1, the 25 
term Q in that case being equal to Q-du 0 -du 3 +3/2 (du 6 -du 4 ) 
. This second form has the advantage of improve discrimi- 
nation between late combustions and combustion misses of 
the engine. 

The term q is still provided by the relation q»6u«(du 0 - 30 
du^+duj-du.^), however, the weighting constants a and p 
used for the calculation of Dc are then different, a remaining 
equal to ji/2V3 and p assuming a value which is determined 
experimentally. 

As indicated in the preamble of the present description, 35 
the formula used for the calculation of the mean gas torque, 
that is Cg-A-D/T^+B/T 2 , is only valid for engines in low- 
level operating conditions. In the case of an engine in 
high-level operating conditions, a correction must be applied 
to the coefficient A in question, as indicated in the patent 40 
referred to in this preamble. For this purpose, the coefficient 
A is replaced by A c =A(l-z) where Zo(fi/fr) 2 . In this equation 
fi =1/T, the instantaneous frequency of the combustions and 
fr is the torsional resonance frequency of the shaft/inertial 
flywheel connection. To apply this correction in the context 45 
of the present invention, the memory 56 will contain the 
constant term fr and the calculation stage 54 will calculate 
fi-l/T and then z=(fi/fr) 2 and then A c =A(l-z), and then, 
finally, the desired value Cg c -A c -D/T 3 +B/T 2 , which, in this 
case, is both valid for low- and high-level operating condi- 50 
tions of the engine and correct for low-frequency perturba- 
tions. 

The invention is not limited to the embodiment described 
with reference to FIG. 1 and to lines I and II of FIG. 3. 
Indeed, it is possible to replace the eight weighted times of 55 
the sets-references of line II of FIG. 3 by the eight other 
weighted times of line III of the same figure. In this case, the 
times which are not taken into account in the calculation of 
Q' c , that is, du x and du 4 , concern six teeth and the other 
times, seven teeth. Consequently, the times duo and du s each 60 
concern two teeth within the angular range separating two 
consecutive high dead centers and five teeth outside of this 
range. The relation connecting Q' c to the times du_ a , du 0 , . 
. . , du 5 , du +1 , represented in line III of FIG. 3, is the same 
as that connecting Q c to these same times represented in line 65 
II of this figure. These two relations are illustrated with the 
signs and weightings in these lines II and III. Similarly, in 


483 
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line IV of FIG. 3, a continuous series is represented, of four 
times of sets-references du 0 , . . . , du 3 each concerning eight 
teeth, the times du 0 and du s each including one tooth outside 
the angular range separating two consecutive high dead 
centers. 

In the case of line III of FIG. 3, the number of sets teeth 
taken into account and those sets not taken into account, as 
well as the weightings and the signs of the times of sets taken 
into account are determined so that the simplified calculation 
term Q'=(du 0 -du 2 -du3+du 5 ) presents a nonzero mean value 
during one measurement horizon and a nonzero frequency 
response at the frequency of this horizon. A constant weight- 
ing coefficient p 1 is then applied to the term Q' so that the 
frequency response of the term D=p' Q* at the frequency of 
the measurement horizon is identical to that supplied for this 
term D obtained by the initial calculation procedure. The 
three parameters indicated above (numbers of teeth, weight- 
ings and signs of the sets) are also determined so that the 
frequency response of the term D=p' Q' at frequencies 
greater than the measurement horizon frequency is as low as 
possible. Moreover, for the calculation of the term q' for 
correction of the low-frequency perturbations, the expres- 
sion defining q' is identical to that of q indicated above, and 
the numbers of the teeth of the sets whose times are 
referenced by du_j and du +1 are the same as those of the 
times referenced du 0 and du s , as appears in line III of FIG. 
3. 

All the above considerations, which concern Q' and p', 
apply to the entire embodiment presented in line IV of FIG. 
3, the value of Q' then being Q'=(du 0 -du 1 -du 2 +du3) 

In the two cases illustrated in lines III and IV of FIG. 3, 
the preparation of the necessary magnitudes for the calcu- 
lation of the angular period T of the combustions of the 
engine will be modified to take into account the fact that the 
measurement horizon extends beyond both sides of the 
angular range concerned. This preparation will consist in 
calculating in addition to the times necessary for the calcu- 
lation of Q' and Q' c , times du 0 , and du s , (for the example of 
fine III) or du 3 , (for the example of line IV) which are times 
du^ du 5 and du 3 which have been truncated and limited only 
to those teeth within this angular range. For this purpose, the 
stage 36 for the measurement of the times of the sets and for 
the identification of their ranks described above will calcu- 
late in addition these truncated times and assign an adequate 
identification to them. The stage 42 for the calculation of the 
period T will be programmed accordingly. 

In a variant, it is noted that the two sets of times du 2 and 
du 3 in the case of lines I, II and III of FIG. 3, and those of 
times du 2 and du^ in the case of line IV, can be grouped into 
a single set by the stage 36. The calculations of T and of Q c , 
which follow, will be modified accordingly. 

It is noted that the angular range defining the width of the 
sets of rank 0 and m' determines the precision of the 
measurement of SCg. When this range is relatively small, 
this precision decreases. When it becomes relatively large, 
the measurement becomes sensitive to the pressures in the 
cylinders. This is due to the fact that the lever arms of the 
connecting rods are no longer negligible and, in this case, the 
set of rank m' of the measurement horizon analyzed before 
(its time is reference du_j) and that of the rank 0 of the 
following horizon (its time is reference du +1 ) are at a 
significant distance from the two consecutive high dead 
centers which determine the theoretical boundaries of the 
measurement horizon during an analysis. To correct this 
negative effect of the broadening of the sets of rank 0 and m\ 
it is possible to introduce a slight angular shift of the 
measurement horizon during an analysis by slightly advanc- 
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ing beyond the boundaries of this measurement horizon with The measurement and calculation stage 30 comprises 

respect to the two high dead centers which define them counting circuits which receive high-frequency timing 

theoretically. The purpose of this is to gain some distance pulses (10 MHz, for example) produced by a quartz clock 

from the effects of the combustion. 32, and it delivers numerical values at the output which are 

It is recalled that the calculation of the negative quality 5 representative of the number of clock pulses counted 

factor 6Cg of the torque Cg can be carried out regardless of between two low-high transitions of the signals with steep 

how the term Cg is calculated, in the simplified or unsim- edges produced b y the shaping stage 26. As far as the 

plified form, and that in addition, it can be made available, pro06Ssirjg of eacb signa i 2 8 produced by one of the long 

as is to the user. For this purpose, an additional calculation ^j. ^ 18 _ 2Q fc concerned in the context of the 

stage 55 will be provided which receives a from the memory describe e d { yalue ^ be b lWQ (flne _ bit 

52, q from the memory 48, A from the memory 56 and T ^ * ^ 

from the memory 44, and which calculates 6Cg=aq-A/T J . . J . . ... . - ft , F 

A circuit 10 is represented in FIG. 4 for the measurement «« virtual tee ^ ' 10 ^ manner the stage 30 for the 

of the corrected mean gas torque in a first step of the low- measurement and the calculation of the periods d, addresses 

and high-frequency perturbations produced by each com- to a buffer memor y 34 successive series of thirty numerical 

bustion of the gaseous mixture in a four-stroke, four- 15 vaJues d * respectively associated with thirty consecutive 

cylinder combustion engine operating at low and high numbers ranging from 0 to 29 defining the rank i of each of 

operating conditions (Cg cl ) and, in a second step, of the the real or virtual teeth within a measurement horizon, 

perturbations affecting the measurement because of the The buffer memory 34 is connected to a calculation step 

particular conditions of the combustions in the engine 36 which is adapted for regrouping, in sets of five, the thirty 

20 instantaneous times d ( . between two reference teeth. In this 

For this purpose, a toothed measurement ring gear 12 is manner, the stage 36 successively produces six passage 

mounted, with integral connection, on the inertial flywheel times of five teeth-reference sets, respectively numbered 

of an engine with electronic ignition. For example, the ring diio, . . . , du 5 as they are produced, the angular value of each 

gear 12 comprises, on its periphery, fifty-six identical, set being x/6, these times and numbers are then transmitted 

regularly spaced, teeth, such as those formed by the full part 25 to a buffer memory 38 placed under the control of a transfer 

14 and the hollow part 16, distributed into two series of control stage 40 which is also adapted to receive the times 

twenty-eight teeth separated by two diametrically opposite du,. 

reference teeth, such as the tooth formed by the solid part 18 The buffer memory 38 is adapted for containing three 

and the hollow part 20, which have a width that is twice of successive series of six times of sets-references: a central 

that of the other teeth. Indeed, the ring gear 12 comprises 30 series which defines the measurement horizon during an 

(2x28+2x2)-60 equidistant references consisting of real or analysis, a front series which precedes this central series in 

virtual teeth of the same module. The two diametrically time and a rear series which follows it. At the arrival of a 

opposite broad teeth serve as reference or origin of index- new time duo, the transfer control stage 40 is adapted for 

ation to allow the numbering of each of the teeth and transferring to the following stages the six times referenced 

especially to identify the tooth b^, which will be defined 35 di^ du 5 of the central series as well as, under the 

below. reference du_ 1} the time du 5 of the last set of the front series 

With the four-strokes, four-cylinder combustion engine and, under the reference du +1 , the time du 0 of the first set of 

mentioned above, it is immediately noted that the angular the series which has just arrived. 

period of the combustions concerns thirty teeth and is equal A first calculation stage 42, which receives these eight 

to half the period of rotation of the crankshaft. 40 successive times, is programmed to select the six times 

A fixed sensor 22 is associated with the ring gear 12, for referenced du 0 , . . . du 5 and to calculate their sum, which is 

example, one with variable reluctance, which is adapted for the instantaneous period T of an angular range of combus- 

the delivery of an alternating signal 24 with a frequency tions. Each new value T so calculated is addressed to a buffer 

proportional to the rate of passage of the teeth of the ring memory 44 where it replaces the value calculated previ- 

gear, that is proportional to the instantaneous speed of the 45 ously. 

flywheel. A second calculation stage 46 which receives these eight 

The angular position of this sensor 22 with respect to the successive times (du_ la du^, . . . , du 3 , du +1 ) is programmed 

indexing teeth 18 at the time when the piston of a cylinder to select the six following times: du. lt du 0 , duj, du 3 , du 5 and 

is at top dead center is known or determined. This allows the du +1 and to combine them according to the relation: 

identification of the tooth do as the one which passes before 50 „ „. . „. . . 
the sensor during the passage of the piston or the given 

cylinder at its high dead center of combustion. The signal The term Q c so calculated is addressed to a buffer memory 

delivered by the sensor 22 is applied to the input of a shaping 48 where it remains until a new value of Q c arrives. The 

circuit 26 which is adapted for the delivery of signals 28 buffer memory 48 is connected to a third calculation stage 50 

with steep edges, having a time equal to the instantaneous 55 which also receives from a memory 52, the weighting 

period d f of the incident signals 24, the index i varying from constant p«n/2V3. The calculation stage 50 is adapted for the 

0 to 29 as the teeth pass before the sensor. Each period d, of calculation of the product D c =pQ c and for its application to 

the signal so produced corresponds to the passage time of a a fourth calculation stage 54. 

tooth, either a solid part or a hollow part, before the sensor This calculation stage 54 also receives two constants A 

22. As far as the incident signals produced by the reference 60 and B from a memory 56, which are experimentally deter- 

teeth 18 are concerned, the shaping circuit 26 transforms mined and the term T from the buffer memory 44. From 

them in the same manner into a signal with steep edges, these two magnitudes, the stage 54 calculates the gas torque, 

having a time that is exactly twice that of the signals corrected for low-frequency pertuibations, according to the 

assigned to the other teeth. The signals 28 are applied to a equation Cg c =AD c /T 3 +B/T 2 . The torque Cg,. so obtained is 

stage 30 for measurement and calculation of the instanta- 65 corrected for low-frequency perturbations, 

neous periods d,- of passage of the real and virtual teeth of the As mentioned in the preamble of the present description, 

measurement ring gear 12 before the sensor 22. the formula used for the calculation of the mean gas torque, 
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that is Cg=A D/T^+B/T 2 is only valid for engines operating 
at low-level condition. In the case of an engine operating at 
high-level conditions, a correction must be applied to the 
coefficient A in question, as indicated in U.S. Pat. No. 
9,111,273 referred to in the preamble. For this purpose, the 
coefficient A is replaced by A^-Afl-z) with z=<fi/fr) 2 . In this 
relation fi*»l/T, the instantaneous frequency of the 
combustions, and fr is the torsional resonance frequency of 
the shaft/inertial flywheel connection. To use this correction 


According to FIG. 5, a mass flow sensor 70 is placed in 
the air-intake manifold. This sensor 70 will be, for example, 
of the type with hot wire or with embedded elastic blade at 
one end. The sensitive elements of these sensors are 
mounted in a bridge, so as to produce first an analog signal, 
which is then transformed into a digital value incorporated 
analog/digital converter. This digital value is then applied to 
a buffer memory 72, which is adapted for holding at any time 


three data elements D„., D„._, and D„ 


which are 


in the context of the present invention, the memory 56 will 10 representative, respectively, of the air flow rates in the intake 

contain the constant term fr and the calculation stage 54 will manifold, during the combustion interval of rank n* during 

calculate fi»l/T and then z=(fi/fr) 2 and then A^oACl-z), and the course of an analysis, and of the two combustion 

then, finally, the desired value Cg c =A c D/T 3 +B/T 2 , which, in intervals which precede it. The data stored in the buffer 

that case, will be valid both for high- and low-level oper- memory 72 is applied to a calculation stage 74 which also 

ating conditions of the engine and corrected for low- is receives from the buffer memory 44, the period T of the 

frequency perturbations. combustion interval during an analysis (this period is prac- 

According to the present invention, the embodiment rep- tically constant over the three successive periods 

resented in FIG. 4 comprises, in addition, a sensor 58 which concerned). The calculation stage 74 is adapted for the 

senses the atmospheric pressure Pa, a pressure sensor 60 Pc successive calculation of the terms (D„._j T) and (D„._ 2 T) 

in the air-intake manifold, and a buffer memory 62 to which 20 and for applying them to a buffer memory 76. The two latter 

are applied the data produced by the pressure sensors 58-60. terms are respectively M„ V1 and M„._ 2 , as already indicated 

These two pressure sensors will be, for example, of the type above (recalling that T=V$N) . Moreover, the data element 

with strain gauge diffused in a silicon chip. First, they D„,_ 2 stored in the buffer memory 72 is applied to a 

produce an analog magnitude which is representative of the calculation stage 78 adapted for calculating (D„._ 2 -T-M 0 ) 

pressure measured. This magnitude is then converted into a 25 from D„._ 2 and from M 0 , the air mass introduced into a 

numerical value by an incorporated analog/digital converter cylinder for a fully opening butterfly valve, which is a 

before being applied to the buffer memory 62. This buffer measured constant, stored in a permanent memory 80. The 

memory 62 is adapted for holding at any time three succes- term (D„._ 2 T-Mp) so calculated is Cp ( „._ 2) . This term is then 
sive measurements of the pressure Pc made during the 

course of three successive combustion intervals. These mea- 30 


surements are the last measurement received Pc ( „^ and the 


applied to a buffer memory 82. 

The constant terms 8, a and h defined above are stored in 
a permanent memory 84. The terras (D„._ 1 -T) and (D„._ 2 -T) 
stored in the buffer memory 76, the term (D„._ 2 *T-M 0 ) 
stored in the buffer memory 82 and the constants 8 and a 
stored in the permanent memory 84 are applied to the 


two preceding measurements Pc^-.j) and Pc^,._ 2) , each 
measurement being identified by its rank n', (n'-l) or (n'-2). 
The buffer memory 62 is placed under the control of a 

calculation stage 64 which communicates the four pressure 35 calculation step 86 which is adapted for the derivation of the 
data which it contains to it. The calculation stage 64 is above-mentioned term H„. In this stage 86, the relation 
adapted for the calculation of the term expressing H„ is: 


40 


from three of the four above-mentioned pressure data 
elements and three constants 8, R„ and a stored in a 
permanent memory 66. The values of these constants 
have been indicated above. 
A last calculation stage 68 receives, on the one hand, the 45 
value Cg,.; produced by the calculation stage 54 and, on the 
other hand, the value of the term H„. calculated by the 
calculation stage 64 as well as a constant h which is also 
stored in the memory 66. The calculation stage 68 is adapted 
for the calculation of the desired value Cg c „ r2 according to 5Q 
the relation 

As indicated above, the magnitude Cg„. 2 thus obtained is 
the completely corrected numerical value of the torque of a 55 
four-stroke, four cylinder internal combustion engine of 
injection, which operates at low- and high-level operating 
conditions, this value being in fact also corrected for low- 
frequency perturbations which affect the engine as well as 
for errors caused by the particular conditions of the com- 60 
bustions in the engine (rate of recirculation of the exhaust 
gases, rate of residual gases, excess of fuel and transitory 
regimens) . As indicated above, the invention can be imple- 
mented using a measurement of the air flow rate of the 
collector, instead of a pressure measurement. The particular 65 
means of production of this variant are represented in FIG. 
5. 


fl>- 6(£>„.. l T+O^_ 2 r)+a(D f ,_ 2 r-Af 0 ) 

This relation is equivalent to the one used above with data 
representative of the pressures in the air collector and of the 
atmospheric pressure. 

During the period of analysis of rank (n*-2), data is 
produced, Pc (n ._ 2 j and D ( „._ 2) , which allow, on the one hand, 
the calculation of the pumping torque Cp (n ._ 2) of the cylinder 
which will be in release during the period of analysis of rank 
n' and, on the other hand, the air mass M (n ._ 2) which will 
participate in the combustion effected during this period of 
rank n'. During the period of analysis of rank (n-1) , data 
P c (n--i) or are produced which allow the calculation 

of the compression torque of the cylinder during 
compression, during the period of analysis of rank n', which 
will be in release during the period of analysis of rank (n'+l); 
this compression torque is proportional to M ( „,_ 1) . 

This pumping torque and this compression torque, as well 
as this air mass participating in the combustion during the 
period of analysis of rank n f are directly dependent on the 
particular conditions of the combustions. Taking them into 
account allows the correction for the effects on the Cg of the 
rate of recirculation of the exhaust gases, of the rate of 
residual gases burned, of the excess fuel and of the transitory 
operating conditions of the engine. 

The invention is not limited to the described examples. As 
a variant, it is possible to calculate Cg from the general 
relation defining D, instead of from the simplified relation 
used in the described example. Similarly, it will be possible, 
in some cases, not to use the corrections pertaining to the 
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low- frequency perturbations and/or the corrections concern- 
ing the high-level engine operating conditions. Moreover, to D _ fa^ 1 ) dj t cos ^. 
take into account the limited variations in the filling coef- 0 
ficienl of the cylinders as a function of the speed of rotation 

of the engine, a correction circuit receiving the mean value 5 2. Method according to claim 1, wherein the number of 

R of this coefficient and the period of the combustions T can reference sets within a measurement horizon is an even 

be provided, which will produce from determined number, and wherein an angular length of these sets are 

experimentally mapping. Such a circuit will be placed equal, and wherein the number of passage times of reference 

between the memory 66 and the calculation stage 64. sets combined to produce Q is a even number equal to at 

We claim- 10 least fouL 

1. Method for the production of a digital value Cg which 3 - Method ^cording to claims 1 or 2, 

is representative of the mean gas torque generated by each wherein the frequency response of the term D-p-Q is 

combustion of the gas mixture in the cylinders of a com- ^ <° or Iess ?™ f ° ur tim f s s f anal y s |f frequency, 

bustion engine, the combustion engine includes measure- 15 A J^f CSSC ^ B ^ Uc ? 1 to h !ha f ™ * I 

ment reference arranged on a ring gear which is integrally . 4 ; Method accordm S «? ^ W f here f in the 

, . • , n u i r .u w • horizon includes a small number of references located on 

connected to an inertial flywheel of the combustion engine bQth gides rf ^ f between tW£) consecutive 

or its crankshaft, means to define at least one reference for losions in the co * bu stion engine. 

indexing the measurement references, a sensor to sense the 5 Method according to claim x wherein the re f eren ce 

passage of the measurement references, which is fixedly 20 mof ameasurement horizon do not have exactly the same 

mounted in a vicinity of the ring gear, said method com- length. 

prising the steps of: 6, Method according to claim 1, wherein the measurement 

producing a primary digital value d, which is represents- horizon * exacd y defined ^ ao 8 ular ran f separating 

tive of an instantaneous passage of time sensed by the 25 »wo explosions m the combustion engine, and wherein the 

sensor of each of the measurement references; ™ mber ° f refer f nce «** » a ° d ' he S1X ^Tn JZ* 

referenced are du 0 , .... du 5 , and the magnitude Q=<du 0 - 

deriving, from the primary digital values d t -, a first sec- di^-dus+dus) and the weighting constant p-n2V3, when the 

ondary digital value T, which is representative of the flxed sets have the same angu i ar length n/6. 

total passage of time before the sensor of each series of 7 Method according to claim 1, wherein the angular 

n references denning an angular range of the combus- 30 rail g e between two consecutive explosions is cut into six sets 

tions in the combustion engine; 0 f identical angular length, and wherein the measurement 

deriving a second secondary digital value D, which is horizon is equal to seven consecutive reference sets du 0 . . . 

representative of a projection onto a phase reference , du 6 , and wherein the magnitude Q is equal to CMUo-dUg* 

fine of the means to define at least one reference for 35 3/2*(du 6 -du 4 ). 

indexing the measurement references, corresponding to « Method for the production of a first numerical value Cg 

an origin of the angular range of the combustions, of an which ■ representative of the mean gas torque generated by 

amplitude of an alternating component of the instanta- each combustion of the gaseous mixture in the cylinders of 

j- r\u * *• a combustion engine and a second numerical value oCg 

neous passage times di of the measurement references a ~. u *" . 6 . e . & 

/, ,, 6 , e r lU u *■ which is representative of the variations of resistant torque 

sensed by the sensor at a frequency of the combustions 40 , c ,• j . ■ . /? 

a^u^u vy iu t wua w applied to a two- or four-cylinder engine, the combustion 

in the combustion engine; cngin< , indudcs measuremcnt rc f ere nces arranged on a ring 

deriving the desired numerical value Cg from the relation gear wmcn j s integrally connected to an inertial flywheel of 

Cg-AD/T^+B/T 2 in which A and B are experimentally me comD ustion engine or its crankshaft, means to define at 

determined constants; least one reference for indexing the measurement references, 

establishing a measurement horizon equal to at least the 45 a sensor to sense the passage of the measurement references, 

angular range between two consecutive explosions in mounted in the vicinity of the ring gear, said method 

the combustion engine; comprising the steps of: 

regrouping the instantaneous passage times di of the producing a primary numerical value d, which is repre- 

references within this measurement horizon in a rela- 50 * ntative of aD instantaneous passage of time sensed by 

tively small number of passage times du, du M of lhe sen *> r of each ° f the references; 

reference sets and of identifying these times du by the denvmg, from the primary numerical values d„ a first 
rank O ... m in each measurement horizon during the secondary numerical value T which is representative of 
f . the total passage of time before the sensor of each series 
course o an an ysis, of ^ refcrences defining an angular range of the corn- 
combining together by addition and subtraction, possibly bustions in the combustion engine; 
weighted, a given number of times having determined a sccoad numerical value D, which is 
ranks, so as to produce a magnitude Q having a mean representative of a projection on a phase reference line 
value which is equal to zero and a frequency response of the mejms tQ defme at least OQe reference for 
which is nonzero at the analysis frequency of the 6Q indexing the measuremcnt references, corresponding to 
combustion engine explosions; and ao or ig m 0 f the angular range of the combustions, of an 
determining the weighting constant p such that the term amplitude of an alternating component of the instanta- 
D-p-Q obtained by a simplified calculation has a fre- neous passage times di the measurement references 
quency response, at least for said frequency of analysis sensed by the sensor at a frequency of the combustions 
and for the combination of the times used to express Q, 65 in the combustion engine; 

which is essentially identical to that obtained by the deriving the desired numerical value Cg from the relation 

term: Cg-AD/T^+B/T 2 in which A and B are experimentally 
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determined constants; and wherein the method for the 
production of said second value oCg comprises the 
sub-steps of 

establishing a measurement horizon equal to at least the 
angular range between two consecutive explosions in 5 
the combustion engine, 

regrouping the instantaneous passage times di of the 
references within this measurement horizon in a rela- 
tively small number of passage times duo, . . . , du m . of 
reference sets and in identifying these times du by the 10 
rank 0 ... m' in each measurement horizon during the 
course of an analysis, the angular values of the time sets 
of rank 0 and m' being equal, 

assigning the reference du_j to the time of the last J5 
reference set of the last measurement horizon analyzed 
earlier and the reference du +1 to the time of the first 
reference set of the next measurement horizon 
analyzed, 

combining four of these times du according to the relation 2 o 
q=(du 0 -du_j+du m .-du +1 ), and 

deriving the desired term 8Cg according to the relation 
SCg=aqA/T 3 , where A is a new constant dependent on 
the angular value of the four given reference sets. 

9. Method according to claim 8, wherein the measurement 25 
horizon is exactly defined by the angular range separating 
two explosions in the combustion engine, and wherein the 
number of reference sets in the measurement horizon is six 
and the angular value of each set is ji/6, and the constant 
a»Jt/2V3. 30 

10. Method according to claim 8, wherein the angular 
range of the reference sets of times reference times di^ and 
du m , is relatively large to produce a more precise term oCg, 
the measurement horizons present a slight phase lead with 
respect to the intervals of high dead centers of the combus- 35 
tion engine. 

11. Method according to claims 6 or 8, wherein, to 
produce a term Q, corrected for low-frequency perturbations 
generated by variations of a resistant torque applied to the 
combustion engine, this term is calculated according to the 40 
relation: 

12. Method for the production of a numerical value Cg c ^ 
which is representative of the mean gas torque, generated by 
each combustion of the gaseous mixture in the four cylinders 
of a four-stroke combustion engine, and corrected for the 
influences of the particular conditions of these combustions, 
said combustion engine includes measurement references 5q 
arranged on a ring gear which is integrally connected to a 
flywheel of the combustion engine or its crankshaft, means 

to define at least one reference for indexing the measurement 
references, a sensor to sense the passage of the measurement 
references, mounted in a fixed manner in a vicinity of the ^ 
ring gear, said method comprising the steps of: 
producing a primary numerical value di which is repre- 
sentative of an instantaneous passage of time sensed by 
the sensor of each of the measurement references; 
deriving, from the primary numerical values d„ a first go 
secondary numerical value T which is representative of 
the total passage of time before the sensor of each series 
of m references defining an angular range of the 
combustions in the combustion engine; 
deriving a second secondary numerical value D, which is 65 
representative of a projection onto a phase reference 
line of the means to define at least one reference for 
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indexing the measurement references, corresponding to 
an origin of the angular range of the combustions, of an 
amplitude of an alternating component of said instan- 
taneous times d ( at the frequency of the combustions in 
the combustion engine; 

deriving an intermediate numerical value Cg from the 
relation Cg=A D/T 3 + B/T 2 in which A and B are 
experimentally determined constants; 

measuring, especially in the air-intake manifold of the 
combustion engine, physical parameters which allow a 
determination of an air mass M pumped into a cylinder 
during the period of combustion of another cylinder, as 
well as a pumping torque Cp required for this purpose; 

storing in a memory the values of these physical param- 
eters or those of the magnitudes M and Cp, during the 
following two periods of combustion; 

combining the data so stored in the memory to produce a 
term H„' for the correction of the torque Cg„. attributed 
to the period of combustion of rank n according to the 
relation H w ,— 5(N„._ 2 +M 2 „._ 1 )+a-Cp (n ._ 2) , where 6 and 
a are constants which depend on the type of combus- 
tion engine; and 

calculating the corrected mean gas torque, generated 
during the period of combustion of rank n, according to 
the equation Cg Crt .=h-Cg M ., where h is a calibration 
constant which is dependent on the type of combustion 
engine. 

13. Method according to claim 12, wherein determination 
of the air mass M includes measuring the pressure Pc in the 
air-intake manifold of the combustion engine and in calcu- 
lating M according to the relation M=Pc R^, where is the 
coefficient, which is essentially constant as a function of the 
speed N of the combustion engine, of filling of the cylinders 
of the combustion engine, and wherein determination of the 
pumping torque Cp includes measuring the atmospheric 
pressure Pa and in calculating Cp according to the relation 
Cp=(Pc-Pa). 

14. Method according to claim 12, wherein determination 
of the air mass M includes measuring the mass flow Da of 
the air in the manifold and in calculating M according to the 
relation M=DaT, where T is the period of the combustions, 
and wherein determination of the pumping torque Cp con- 
sists in calculating Cp according to the relation Cp=(Da T- 
M D ), where M 0 is the air mass introduced into a cylinder by 
a fully opened butterfly valve. 

15. Method according to one of claims 12, 13, or 14, 
wherein the term D of the relation defining Cg is obtained by 
the performance of the following steps: 

establishing a measurement horizon equal to at least the 
angular range between two consecutive explosions in 
the combustion engine and dividing this measurement 
horizon in a relatively small number of reference sets; 

calculating the passage times of each of these reference 
sets before the sensor and ranking them in each mea- 
surement horizon; 

combining together, by addition and subtraction, a given 
number of times having predetermined ranks, so as to 
produce a magnitude Q having a zero mean value and 
a nonzero frequency response at the analysis frequency 
of the explosions of the motor; and 

determining a weighting constant p such that the term 
D-p-Q obtained by a simplified calculation has a fre- 
quency of response at least for said analysis frequency 
and for the combination of the times used to express Q, 
which is essentially identical to that obtained by a term 
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rf r ooB(i'-2Ji/m). 
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express Q having been determined so that the term D so 
obtained by simplified calculation has a frequency 
response, at least for said analysis frequency, which is 
identical to that obtained for the term: 


16. Method according to claim 15, wherein the angular 
value of each reference set is n/6. 

17. Method according to claim 15, further comprising a 
term D c -p Q,. corrected for the low-frequency perturbations 
which affect the combustion engine, each reference set 

measures n/6, the constant p=Ji/2V3 and the value of the 10 representative of the mean gas torque generated by each 


19. Device to produce a first numerical value Cg which is 


term Q so corrected is Qc-(2duo-du. 3 du 2 -du3-du +1 +2du 5 ), 
with du.., and du +1 , the respective passage times of the sets 
arranged at a given instant on both sides of said angular 
range. 

18. Device for the production of a numerical value Cg 
which is representative of the mean gas torque generated by 
each combustion of the gaseous mixture in the cylinders of 
a combustion engine at low-level operating conditions, the 
combustion engine includes measurement references 
arranged on a ring gear which is integrally connected to an 
inertial flywheel or a crankshaft of the combustion engine, 
means to define at least one reference for indexing the 
measurement references, a sensor to sense the passage of the 
measurement references which is mounted in a fixed manner 
in a vicinity of the ring gear, said device comprising: 
calculation means to produce a primary numerical value 
d, which is representative of an instantaneous passage 
of time sensed by the sensor of each of the measure- 
ment references; 
calculation means to derive from the primary numerical 
values d,. a first secondary numerical value T which is 
representative of a total passage time before the sensor 
of each series of n references defining an angular range 
separating two consecutive combustions in the com- 
bustion engine; 
calculation means to derive a second secondary numerical 
value D which is representative of a projection, onto a 
phase reference line of the means to define at least one 
reference for indexing the measurement references, 
corresponding to an origin of the angular range of the 
combustions, of an amplitude of an alternating com- 
ponent of the instantaneous passage times di of the 
measurement references sensed by the sensor at a 
frequency of the combustions in the combustion 
engine; 

calculation means to derive the desired numerical value 
Cg from the relation Cg=A D/T a +B/T 2 in which A and 
B are experimentally determined constants, contained 
in a memory; 

processing and calculation means to establish a measure- 
ment horizon equal to at least said angular range, to 
regroup the instantaneous passage times d, of the 
references within this measurement horizon in a rela- 
tively small number of passage times duo, . . . , du m of 
references sets and to define these times du by their 
rank (0 ... m) in each measurement horizon; 

calculation means to combine together, by addition and 
subtraction, a given number of set times having prede- 
termined ranks, so as to produce a magnitude Q having 60 
a zero mean value zero and a nonzero frequency 
response at the analysis frequency of the explosions of 
the combustion engine; 

calculation means to produce the term D-p-Q from the 
value Q calculated earlier and a weighting coefficient p 
determined experimentally and stored in a memory, the 
value of p and the combination of the times used to 
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combustion of the gas mixture in the cylinders of a com- 
bustion engine operating at low-level condition, and a sec- 
ond numerical value 5Cg representative of the variations of 
the resistant torque applied to such an engine comprising 
two or four cylinders, said combustion engine includes 
measurement references arranged on a ring gear which is 
integrally connected to a flywheel or a crankshaft of the 
combustion engine means to define at least one reference for 
indexing the measurement references a sensor to sense the 
passage of the measurement references, which is mounted in 
a fixed manner in a vicinity of the ring gear, said device 
comprising: 

calculation means to produce a primary numerical value 
d,- which is representative of an instantaneous passage 
of time sensed by the sensor of each of the measure- 
ment references; 

calculation means to derive from the primary numerical 
values d, a first secondary numerical value T which is 
representative of a total passage time before the sensor 
of each series of n references defining an angular range 
separating two consecutive combustions in the com- 
bustion engine; 

calculation means to derive a second secondary numerical 
value D which is representative of a projection, onto a 
phase reference line of the means to define at least one 
reference for indexing the measurement references 
corresponding to an origin of the angular range of the 
combustions, of an amplitude of an alternating com- 
ponent of the instantaneous passage times d, of the 
measurement references sensed by the sensor at a 
frequency of the combustions in the combustion 
engine; and 

calculation means to derive the first desired numerical 
value Cg from the relation Cg-AD/T'+B/r 2 in which 
A and B are experimentally determined constants con- 
tained in a memory, wherein said device to produce 
said second value 6Cg further comprises, 
processing and calculation means to form a measure- 
ment horizon equal to at least one angular range 
between two consecutive explosions in the combus- 
tion engine, to regroup the instantaneous passage 
times d, of the references within this measurement 
horizon in a relatively small number of passage times 
du 0 , . . . , du m of reference sets and to identify these 
times du by the rank 0 ... m' in each measurement 
horizon during the course of an analysis, the angular 
values of the reference sets of the times of rank 0 and 
m' being equal, 
means to store in memory the times du 0 , . . . , du^' as 
well as, on the one hand, the time of last reference set 
of the last measurement horizon analyzed previously 
and to assign it the reference du_j and, on the other 
hand, the time of the first reference set of the next 
measurement horizon analyzed and to assign it the 
reference du +1 , 
means to produce a term q according to the relation 
q=(du 0 -du_ 1 +du m .-du +1 ), and 
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means to derive the desired term 6Cg according to the 
relation dCg-aq-A/T 3 where A is a constant which 
depends on the angular value of the four reference 
sets concerned. 
20. Device for the production of a numerical value Cg 5 
which is representative of the mean gas torque generated by 
each combustion of the gaseous mixture in the four cylinders 
of a four-stroke combustion engine at low-level operating 
conditions, said value being corrected for influences of the 
particular conditions of these combustions, said combustion id 
engine includes measurement references arranged on a ring 
gear which is integrally connected to an inertial flywheel or 
a crankshaft of the combustion engine, means to define at 
least one reference for indexing the measurement references, 
a sensor to sense the passage of the measurement references is 
mounted in a fixed manner in a vicinity of the ring gear, said 
device comprising: 

calculation means to produce a primary numerical value 
d,- which is representative of an instantaneous time of 
passage time sensed by the sensor of each of the 20 
measurement references; 
calculation means to derive from the primary numerical 
values d,- a first secondary numerical value T which is 
representative of a total passage time before the sensor 
of each series of m references denning an angular range 
separating two consecutive combustions in the com- 
bustion engine; 
calculation means to derive a second secondary numerical 
value D which is representative of a projection, onto a 
phase reference line of the means to define at least one 
reference for indexing the measurement references, 
corresponding to an origin of the angular range of 
combustions, of an amplitude of an alternating com- 
ponent of the instantaneous passage times d,- of the 
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measurement references sensed by the sensor at a 
frequency of the combustions in the combustion 
engine; 

calculation means to derive an intermediate numerical 
value Cg from the relation Cg-A-D/T 3 +B/T J , where A 
and B are experimentally determined constants, con- 
tained in a memory; 

means including sensors installed in an air-intake mani- 
fold to determine an air mass M pumped into a cylinder 
during the course of the period of combustion of 
another cylinder; 

means to determine the pumping torque Cp required for 
this purpose; 

means to store in memory the data M and Cp or the 
constituent elements of these data during the two 
following combustion periods; 

means to combine this data or its constituent elements to 
produce a term H. of correction of a torque Cg„. 
pertaining to the period of combustion of rank n 
according to the relation, 

H„'— S-(M n ._2+M 2 „,_ 1 )+a.-Cp in ._2 ) , in which 8 and a are 
constants which are dependent on the type of combus- 
tion engine, 

calculation means to calculate the corrected mean gas 
torque, generated during the course of the combustion 
period of rank n', according to the relation, 

Cg^hCg^+H^., in which h is a calibration constant 
which is dependent on the type of combustion engine; 
and 

the constants 5, a and h being stored in a permanent 
memory. 
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